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Abstract 

In this paper we discuss a quantity that can be derived from the hadronic spectrum 
- the ratio between the color hyperfine splitting of a K meson and that of a £ 
baryon. It is shown that within the constituent quark model this ratio depends only 
on the ratio of contact probabilities in the hadrons. We compute this ratio assuming 
several phenomenological potential models, and show that the best agreement with 
data is obtained from the Cornell potential - Coulomb + linear. Thus the analysis 
of color hyperfine interactions corroborates the preference for the Cornell potential 
based on quarkonium spectra. 

Key words: QCD, Quark potential models, Color hyperfine interaction 



1 Introduction 



The recent measurement of the E& baryon masses |T] provides us with a rare 
opportunity to better our understanding of the interaction between a heavy 
quark and a light diquark at low energies. In this paper we study the dynamics 
of such a baryon by comparing its properties with those of a meson which has 
similar quark content. We focus on their color hyperfine splitting (HF), which 
is defined as the mass difference between 2 similar hadrons that differ by one 
unit of spin (e.g. Ak = Mk* — Mk)- According to the available experimental 
data, the ratio between the HF splitting of a uq meson and that of a udq 
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baryon (where q stands for s, c or b) is almost independent of the heavy quark 
flavor: 

A K /A E « A D /A Sc « A B /A Sb « 2.1 (1) 



This quantity was first discussed in [2], but so far we have not been able to 
explain this result from first principles. However, when using phenomenolog- 
ical models this ratio takes on a new meaning - it will be shown below that 
within the constituent quark model, this ratio is equal to the ratio of contact 
probabilities between quarks inside the hadron. For example: 

_A (jj\5(f u -f s )\^} K 



The interesting feature of this quantity is that it is determined by two param- 
eters only: the shape of the potential that binds the quarks together, and the 
quark mass ratio. While the constituent quark masses can be extracted with 
sufficient accuracy from relations between other hadron masses j3], the shape 
of the potential is still under continuing research. This subject was studied in 
the past mainly through the analysis of the quarkonium spectrum and the re- 
sults favored a Coulomb + linear potential [8|9lfT0] . but this result still awaits 
a solid theoretical justification. 

We therefore suggest to use this ratio as a new means to test the possible 
quark confinement models. We computed the ratio of contact probabilities 
assuming various phenomenological models: harmonic oscillator, Coulomb in- 
teraction, linear potential, logarithmic potential and the Cornell potential (lin- 
ear +Coulomb), and compared the different results with the available data for 
the HF splitting ratio. 



2 The constituent quark model 



Low-energy constituent-quark models have been surprisingly successful in pre- 
dicting hadron masses [3f4"] . In such models the hadrons are considered as 
bound states of constituent quarks in a confining potential, and their total 
mass is given by the Sakharov-Zeldovich formula [5|6f7] : 

M hadron = Y,m l + J2V^ F , (3) 

i i<j 



where m ; is the constituent quark mass and V^ F is the color hyperfine (HF) 
interaction between the quarks i and j. The HF interaction is described as 



2 



the product of the magnetic moments of the quarks, in analogy with the 
electromagnetic hyperfine interaction, and is given by the expression: 



i>j 



m,:m 



(4) 
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where \ are the SU(3) generators, rrii and o~i are the mass and spin of the z'th 
quark, is the distance between the quarks i and j, and t>o is the coupling 
constant (e.g., in the case of one-gluon exchange between quarks and anti- 
quarks we have v X q ■ Xg = — g~~~)- 

We note that this is a contact interaction and that the formula is correct for 
S-wave hadrons only. We also assume that the same formula and parameteri- 
zation apply for both mesons and baryons. 



2.1 Extracting constituent quark masses 



We begin by demonstrating calculations with which constituent quark masses 
can be extracted from the hadron spectrum. The first example is the extraction 
of quark mass differences from A baryon masses: 

A c - A = m c - m s , (5) 

where we used the fact that the HF interaction between the heavy quark and 
the light quarks must be zero, due to the anti-symmetry of the A baryons 
under exchange of the light quarks. 

We can also use the expression for the meson HF splittings 

M K * - M K = Av ^^(5(r u ,)) (6) 
m u m s 

to extract the ratio between the s-quark and c-quark masses: 



M K * - M K \. X- m c ' 



(7) 



2 We use here the normalization under which the spin products of a vector meson 
and a pseudoscalar are {(f u ■ as)s=i = 1 and (cT u ■ (Js) s =o = —3 
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where we implicitly assumed that the interaction strength vq is equal in both 
cases and that the contact term does not vary greatly. 

Given the simple form of this model, and the fact that we simply ignored the 
dependence on the contact probabilities, the consistency of the predictions 
given by this approach is surprising. The following values for the constituent 
quark masses can be extracted from these calculations [3]: 



m c 



360 MeV 
540 MeV 
1710 MeV 
5050 MeV 



2.2 The ratio between HF splitting of baryons and mesons 



We now turn to quantities that involve both mesons and baryons, where the 
contact probabilities are no longer expected to be similar. We choose to an- 
alyze the ratio between two HF splitting values (Mr* — M^)/(Ms* — Ms), 
because as will be seen below, some cancellations occur that simplify the final 
expression. 

The K Meson HF splitting is given by (JH]). The calculation of the HF splitting 

of the £ baryon is easy because it gets no contribution from the ud diquark. 

3 

The explanation for that is as follows: The total spin of the £* is -, so all 

the quark pairs must be at relative spin 1. The £ has isospin-1, therefore the 
2 light quarks must also have relative spin 1, and the contribution of the ud 
diquark to the HF splitting is zero. 

We are left only with the interaction between the light quarks and the heavy 
one. After some spin products calculations we get the following expression for 
the £ HF splitting: 



m„rn s 



X u -K(6(r us )) k (9) 

J baryon 



Combining equations ([HD and (TTU|) . and the relation: 

[Aj • ^j\meson 2[Aj • ^j\baryon (10) 
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we reach the result mentioned above: 

meson /, , \ 

A if/ A 2 - o 77TT7 TfT l 11 ; 

3 (^|<H r «J I VV baryon 



The advantage of this expression over relations such as ([7]) (besides the fact 
that we don't assume anything about the wave function), is that here we have 
the same quark content in all the interactions involved, and it seems safer 
to assume that interaction parameters (confinement strength, HF interaction 
coupling constant) are identical and therefore can cancel out in the ratio. 

Similar expression are obtained when the s quark is replaced with c or b quarks. 
The available experimental data is given in table [TJ 



Meson 


Baryon 


7713/7711 


HF splitting 
ratio 


K 


£ 


1.33 


2.08 ±0.01 


D 




4.75 


2.18 ±0.08 


B 


s 6 


14 


2.15 ±0.20 



Table 1 

Experimental data of meson / baryon HF splitting ratios. 

An interesting comment is that if we simply assumed that the contact proba- 
bilities are inversely proportional to the number of quarks in the hadron, then 
the expected HF splitting ratio would be 2, which is not very far from the 
measured data. 



3 Confining potentials 



Following is a list of phenomenological potentials that are often used in the 
literature as models for the confining potential: 

• Harmonic oscillator (HO) - For sufficiently small deviations from equilib- 
rium, almost every potential looks like a harmonic potential. Of course, 
as the deviation from equilibrium gets larger, the potential is less reliable. 
Nevertheless, it is useful as a proof of concept and for establishing orders of 
magnitude. An additional advantage is that expectation values for 3-body 
problems can be calculated analytically. 

• Coulomb interaction - The leading order QCD diagram involves a single 
gluon exchange, which gives a Coulomb-like potential. 
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• Linear potential - The asymptotic behavior of the QCD confining potential 
at large distances is linear. 

• Coulomb + Linear (Cornell potential) - The Coulomb + linear potential is 
referred to as the Cornell potential, named after the team which studied it 
in the late 70's [8|9] . This potential is a sum of the two asymptotic limits, 
and is supported by theoretical calculations using Wilson loops [11] and 
renormalon cancellation [12] techniques. 

• Logarithmic potential - The logarithm function bears a resemblance to a 
Coulomb+linear function, and is also used as a phenomenological confining 
potential. Several aspects of the hadronic spectrum are consistent with a 
logarithmic potential model. 

The problem of choosing the most effective confining potential model was stud- 
ied in the past mainly through the analysis of quarkonium spectra and leptonic 
decays (TU], and the results favored the Cornell potential. We addressed this 
issue by computing the contact probability ratio from Eq. ffTTl) given each of 
the above models, and comparing with the available measurements. 

The quark mass ratios used in the following calculations were taken from Eq. 
([S]). Due to the error bars in the experimental data, the level of accuracy 
required from the computation of contact probabilities is in the order of 1%. 

3.1 Harmonic oscillator 

In this model the potential that binds the pairs of quarks takes the form 

^E^VA^Xn-rl) 2 (12) 
i<j z 

where V^ is an effective coupling constant, and the product of the SU(3) 
generators depends on the specific quark configuration in question (meson or 
baryon). 

3.1.1 HF splitting in a harmonic oscillator model - mesons 

The ground-state wave function of a harmonic oscillator for a system of two 
masses mi, rri2 is 

Mr) = (^)le-^ 2 (13) 
where we use the definitions 
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mi + m 2 ' m r 



_ m T ijj r 
a r = 



h 

leading to 

{MS(ru,Mo) = \Mo)\ 2 = (^) 3 (15) 



7T 



3.1.2 HF splitting in a harmonic oscillator model - baryons 

The three body HO Hamiltonian, with mi = m 2 and (Ai • A3) = (A2 • A3), can 
be decoupled by a simple transformation of coordinates: 

f=fi-f 2 -R = ^ - 77(^1 + ^"2) . (16) 



These are the two orthogonal modes described in fig. [TJ 




Fig. 1. The two orthogonal modes of the 2+1 harmonic oscillator system 
We now define the reduced masses and frequencies of these modes: 



_ / / by — 

_ 2mim 3 2 _ Vfo -> - 

mR = 2m7T^ UR = m^ 2(Al • As) 

The new Hamiltonian is a sum of two commuting Hamiltonians; therefore the 
wave function will be a product of the eigenstates of these Hamiltonians: 

MR, r) = <Mr> r0 (r) = (— )i e 4K* 2 +«V) (18) 

7T 



7 



using the definition given above - aj — v /; . 
The calculation of the 5-function expectation value is as follows: 



^1^0(2)11^0(^1 



(^) 3 ( / dxe-<^y 

1 d r a R 



(19) 



We can now plug the masses from ([8]) into the expressions we reached, and 
compute the HF splitting ratios for the meson and baryons under consideration 
(table [2}. 



Meson 


Baryon 


HF splitting 


Deviation 






ratio 


from data 


K 


S 


1.65 


21% 


D 




1.62 


26% 


B 


s 6 


1.59 


26% 



Table 2 

Prediction of meson / baryon HF splitting ratios for the harmonic oscillator model. 



3.2 Coulomb interaction 



The three-body problem with Coulomb interactions has no analytic solution, 
and we have to use variational methods in order to calculate the contact 
probabilities. The ansatz we chose was introduced by Coolidge and James 
[131111] as part of a modification of the so-called Hylleraas method, and it uses 
the isoperimetrical coordinates: 



u = r i2 + r 23 - ri 3 
v = r 12 ~ r 23 + ria 

w = -ru + r 23 + 7*13 (20) 
The integration measure in these coordinates is: 
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7T 2 

^ r i3^ r 23 = -^-{u + v){u + w){w + v)dudvdw 
and the kinetic terms take the form: 

d 2 d 2 d 2 



Pi 



2m\ mi{u + v )(v + w) 



( 2 1 (21) 

^du 2 dw 2 dudw' 

+ v{u+v+w) h + {w - u) ^u~L) 

d ' 

+ (u + 2v + w) — 
ov 

The kinetic terms for the second and third particles can be obtained by a 
permutation of the parameters u,v and w. 

On the basis of the solution for the 2 body problem with Coulomb interac- 
tion, and the requirement for a wave function which is symmetric under the 
exchange of the two light quarks (the E has isospin-1), we choose the following 
ansatz: 

ip(u,v,w) = e -|W«+^)+/3M) x (22) 
C k i m \L k (au)Li(av) + Li(au)L k (av)\L m (f3w) 

k,l,m 

where a , (3 and C k i m are variational parameters, and L k {x) are Laguerre 
polynomials of degree [k, 0] . 

The extraction of the parameters was performed using the Ritz variational 
method. Summing up to polynomials of rank N = k + l + m = 5, we were able 
to reproduce calculations of given in [T3] with a deviation of less than 1%. 

Given the ground state wave function we can calculate the contact probability 



baryon 2 

dv\ip(0,v,0)^ (23) 



leading to the results shown in table El 
3.3 Linear potential 

The eigenstates of a 2-body Hamiltonian with a linear potential 

V lmear {r) = r/a 2 (24) 
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Meson Baryon HF splitting Deviation 
ratio from data 

K S 5.07 ±0.08 144% 

D S c 5.62 ±0.02 158% 

B £& 5.75 ±0.01 167% 

Table 3 

Prediction of meson / baryon HF splitting ratios for the Coulomb interaction model. 
The error bars are obtained by considering a 5% error in the quark mass. 



and reduced mass fx = 1 are 



. , , 1 ..r2r-2a 2 E n 
tfj n {r) = -At - 

r 1 (2a) z/6 



(25) 



where Ai is the Airy function. The energy E n is determined by the boundary 
condition that the Airy function should equal zero as r — > 0: 



with a n the n'th zero of the Airy function. 

Once again, the 3-body problem requires usage of variational methods. In 
order to avoid the difficulty of performing many integrations of Airy functions 
and their derivatives, we tried to use the modified Hylleraas ansatz (Eq. |2"51) 
for the linear potential problem. Testing this ansatz for the 2-body problem, 
with the wavefunction given by 

JV 

Hr) = T, C ^ ar L k (ar) , (27) 
k=i 



we found out that it was enough to take polynomials of degree iV„ 



5 for 



the variational solution to converge to the analytic one (see plot 2(b)) 



Following the success with the 2 body problem we felt confident enough to use 
the same software to calculate the contact probability ratio between mesons 
and baryons (see table HI). 
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(a) 
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CO 
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CO 

o 

• 
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CO 
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0.5 







-0.5 



N=1 
N=3 
N=4 
N=5 








2 

r 



(b) 

Fig. 2. Plot (a) shows the analytic solution for a 2 body ground state with a 
linear potential (r is given in the units where c = h = /j, = a = l). Plot 
(b) shows the accuracy of the solutions reached using the variational method 
(1 - \i '(r) variational? ' /\^{r) ana i ytic \ 2 ) with different polynomial degrees N. 
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Meson Baryon HF splitting Deviation 
ratio from data 



K S 1.88 ±0.06 14% 

D S c 1.88 ±0.08 10% 

B £& 1.86 ±0.09 13% 

Table 4 

Prediction of meson / baryon HF splitting ratios for the linear potential model. The 
error estimates are based on comparison of the N max = 5 results with the analogous 
results for N max = 4. 

3.4 Cornell potential - Coulomb + linear 



Vcornell{r) = -- + ^ (28) 

r cr 



The Coulomb + linear potential combines the behavior of the color interaction 
in the two asymptotic limits, and it is theoretically supported by calculations 
using several techniques [TT][T2"] . This model was used by the Cornell group 
[8] and proved very successful in explaining various aspects of the hadronic 
spectrum. 

Unlike the previous models, in this case the interaction strength does not 
cancel, and the HF splitting ratio depends both on the mass ratio m 3 /mi and 
on an additional parameter: 

K = K(m u a)i (29) 



The parameters extracted from the charmonium spectrum in [9] are equivalent 
to the range 0.25 < K < 0.45. a and k include the product of the SU(3) 
generators; therefore their values are different for mesons and baryons. Using 



Eq. dlDJ) we reach the relations K baryon = \K meson and a baryon = y/2, 



Once again, we tried to solve the 2 body problem using the Coulomb ansatz, 
with satisfactory results - we were able to reproduce calculations for meson 
expectation values given in [91 table I] with deviation less than 1% (N max = 5). 
Using this method to calculate the meson/baryon HF splitting ratio for the 
Cornell potential with 0.2 < K < 0.5 we got the following results: 
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2.04 < A K /A S < 2.28 
2.08 < A D /A Sc < 2.40 
2.07 < A B /A Ei) < 2.43 

The results for the specific value of K = 0.28 are given in Table 

Meson Baryon HF splitting Deviation 
ratio from data 

K S 2.10 ±0.05 1% 

D S c 2.16 ±0.07 1% 

B S fe 2.17 ±0.08 1% 

Table 5 

Prediction of meson / baryon HF splitting ratios for the Coulomb + linear potential 
model with K = 0.28. The error estimates are based on comparison of the N max = 5 
results with the analogous results for N max = 4. 

3.5 Logarithmic potential 

V log (r) = C\og(-) (30) 
r 

One of the consequences of taking the logarithm as the function that describes 
the confining potential is that the energy level spacings are independent of the 
constituent quark masses. This phenomenon is observed when comparing the 
charmonium and the bottomonium spectra, and this is one of the strongest 
reasons why the logarithmic potential is considered a strong candidate for the 
confining potential. 

In order to compute the HF splitting ratio with this potential, we used the 
same technique described above for the Coulomb and linear potentials. The 
theoretical results for the IS" quarkonium state given in [IUJ tables 6-7] were 
accurately reproduced. 

Unfortunately, when using the modified Hylleraas ansatz (T2"3"|) the computation 
of the matrix elements (3)\ \og(r)\ipj(a, (3)) requires careful treatment at 

the limit (a — (3) — > 0, otherwise the integrals diverge. In order to avoid this we 
chose to take into account only a single variational parameter (i.e. we assumed 
a = (3). In order to test whether this modification might generate large errors, 



13 



we checked its effect on the results of the previous models, and the corrections 
did not exceed 1%. The final results for the logarithmic potential are given in 
table O 



Meson 


Baryon 


HF splitting 


Deviation 






ratio 


from data 


K 


E 


2.38 ±0.02 


14% 


D 




2.43 ±0.02 


11% 


B 


s 6 


2.43 ±0.01 


13% 



Table 6 

Prediction of meson / baryon HF splitting ratios for the logarithmic potential model. 
The error estimates are based on comparison of the N max = 5 results with the 
analogous results for N max = 4. 

4 Conclusions 

In this research we tried to explain the measured ratio between the HF split- 
tings of mesons and baryons using the constituent quark model. Thanks to the 
simple fact that the £ baryon is an isospin-1 particle with quark content similar 
to that of a K meson, most of the degrees of freedom in the Sakharov-Zeldovich 
model are eliminated, and the HF splitting ratio was shown to depend only 
on the quarks' contact probabilities. 

Assuming different confining potentials we calculated the expected values for 
this measurement using variational methods. The software that implements 
these calculations reproduced similar calculations found in the literature with 
a deviation less than 1%. 

The results, which are summarized in Fig. [21 show that this quantity enables 
us to distinguish between the different confinement models, and that the ex- 
perimental data favor a Cornell potential, similarly to the results obtained by 
previous experiments. 
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Fig. 3. The K-T, HF splitting ratio as a function of the quark mass ratio assuming 
the different confining potentials. Most of the Coulomb model graph is outside the 
range of this plot. The area marked in gray shows the results that correspond to 
a Cornell potential with K values between 0.2 and 0.5. One can see the good fit 
between the data and the predictions obtained with the Cornell potential (K=0.28). 
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